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Abstract 

The Alcubierre spacetime was simulated by means of a Tamm medium which is asymptotically identi- 
cal to vacuum and has constitutive parameters which are continuous functions of the spatial coordinates. 
Accordingly, the Tamm medium is amenable to physical realization as a nanostructured metamaterial. 
A comprehensive characterization of ray trajectories in the Tamm medium was undertaken, within the 
geometric-optics regime. Propagation directions corresponding to evanescent waves were identified: these 
occur in the region of the Tamm medium which corresponds to the warp bubble of the Alcubierre space- 
time, especially for directions perpendicular to the velocity of the warp bubble at high speeds of that 
bubble. Ray trajectories are acutely sensitive to the magnitude and direction of the warp bubble's ve- 
locity, but rather less sensitive to the thickness of the transition zone between the warp bubble and its 
background. In particular, for rays which travel in the same direction as the warp bubble, the latter acts 
as a focusing lens, most notably at high speeds. 

PACS numbers: 04.40.Nr, 41.20.-q, 42.15.-i 
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1 Introduction 

Metamaterials provide opportunities to study general-relativistic scenarios which would otherwise be either 
impractical or impossible to explore [1]. This may be achieved by exploiting the formal analogy that exists 
between light propagation in vacuum subjected to a gravitational field and light propagation in a certain 
nonhomogcneous bianisotropic medium, known as a Tamm medium [2, 3, 4]. The practical creation of Tamm 
mediums is beginning to look like an increasingly realistic proposition, as rapid developments are made in 
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the science of nanostructured metamaterials. In recent years, metamaterial-based simulations of black holes 
[5], cosmic strings [6], de Sitter spacetime [7, 8] and wormholes [9] have been proposed, for example. 

As a useful Gedankenexperiment, the spacetime of the Alcubierre metric has generated considerable 
interest since its introduction in 1994 [10, 11, 12]. The spacetime is characterized by a warp bubble which 
moves with respect to an asymptotically flat background at a speed v'gCo, where Co is the speed of light in 
the absence of a gravitational field and v'^ > 0. By effectively contracting the spacetime ahead of the warp 
bubble and expanding the spacetime behind it, arbitrary relative speeds v'^ may be achieved, in principle. In 
contrast, within the local environment of the warp bubble, its speed is necessarily subluminal; i.e., v'^ < 1. 

Serious obstacles stand in the way of a physical realization of the Alucbierre spacetime, stemming from 
violations of various energy constraints [13, 14]. Even the subluminal regime is associated with negative 
energy densities [15]. However, negative energy densities are not unprecedented in astrophysics. For exam- 
ples, the emission of Hawking radiation by a black hole is accompanied by a flow of negative energy [16], 
and the construction of wormholes relies on negative energy density [17]. Nor is negative energy density un- 
precedented within the realm of metamaterials: propagation of (monochromatic) plane waves with negative 
phase velocity [18] — which is intimately related to the phenomenon of negative refraction [19] that certain 
metamaterials have been shown to support in experimental observations [20] — also involves negative energy 
density [21], at least in the absence of dissipation [22]. 

In this communication, we present a simulation of the Alcubierre spacetime, in the form a Tamm medium 
which is physically realizable, in principle. Under the geometric-optics approximation, a comprehensive 
characterization of light ray trajectories through the Tamm medium is provided. 

As regards notation, 3-vectors are underlined and unit 3-vectors are additionally distinguished by a caret, 
whereas 3x3 dyadics [23] are double underlined. The identity dyadic is written as /. The scalar constants 
Co and fj,o denote the permittivity and permeability of vacuum in the absence of a gravitational field, and 

Co = 1/ ^Co/Xq. 



2 Tamm medium for Alcubierre spacetime 

The Alcubierre spacetime is specified by the line element [10] 

ds^ = -c^dt'^ + [dx' ~ v'J{r'^)c„dt'f + dy'^ + dz'^, (1) 
for the case where the warp bubble travels at relative velocity j/^ = v'^x . Herein the scalar function 

tanh [cr (r^ + R)] - tanh [a (r^ - R)] 



fir's) 



2 tanh (cri?) 



(2) 



is unit-valued at the origin and decays uniformly to zero as r'^ — >■ oo, where x'g{t') = VgCgt' and the time- 
dependent translated displacement 



r's{t') - ^[x' ^x',{t')f + y'^ + z'' 



(3) 



The scalar parameter i? > is a measure of the warp bubble radius whereas cr > is a measure of the 
thickness of the transition zone between the warp bubble and its background. For cr 3> i?, the function / 
has an approximately top- hat profile which propagates along the positive x' axis at relative speed v'^. 

In order to eliminate the time-dependency which enters via x'^it'), let us introduce the spacetime coordi- 
nates [23] 

t = v it' 



X ~ V [x' — v'^Cat') 

y^y' 

z = z' 



> , 



(4) 
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with the scalar quantity 



This change amounts to a Lorentz transformation [23]. The hue element (1) may then be expressed as 



(5) 



[fir.) - If ivif - 1 1 , {[/(..) - 1] iv'^f - l} 



1 - 

+ {v'sf {firs) [firs) (v'sf 



iv's)' 



dx^ + dy^ + dz^ 



Codtdx 



wherein 



(6) 



(7) 



is independent of t. 

Per the noncovariant approach of Tamm [2, 3, 4], electromagnetic fields in curved spacctime may be 
described by the constitutive relations 



^(r, t) = fc) • KirL. t) - ^^,T (r) x ^(r, t) 
B{r, t) = yi^r (r) X E{r, t) + fiol fc) • H{r, t) 



(8) 



of an equivalent medium in flat spacctime, using SI units, with r — xx + yy + zz. The components of the 
3x3 dyadic 7(r) and the 3- vector T_{r) are defined in indicial notation as 



lim = 



gom 

900 



500 



(£, me {1,2,3}), 



(9) 



where g denotes the determinant of spacctime metric gai3, (a, /? G {0, 1, 2, 3}), prescribing the curved space- 
time. The fictitious Tamm medium represented by the constitutive relations (8) is spatiotemporally local, 
and generally nonhomogeneous and bianisotropic. The sign of the square root term in the definition of ^£m 
is selected such that the metric for vacuous Minkowskian spacctime is represented by the dyadic 7 = /. 

For the case of Alcubierre spacctime characterized by the line element (6), the definitions (9) deliver 



liz) = XX + 



1 - iO' 



7-2 (yy + ll) 



and 



rfc) = 



l-[l-f{rs)r{v',) 
/(nsK{l + [l-/(r.)] jv'f] 



(10) 



(11) 



As we are interested in a physically realizable metamaterial that represents the Alcubierre spacctime, the 
limits 



lim 7(i.) = x£+ 1 — (Wg)^ 
\l\^o = I 

lim r(r) = v' X, 



[yy + zz , 



lim 7(r) = /, 
lim r(r) = 0, 



(12) 
(13) 
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bear considerable promise. Thus, the Tamm medium is like a gravitation-free vacuum for large values of |r| 
whereas its constitutive parameters remain bounded at small values of \r\. 

The nontrivial constitutive parameters included in 7(r) and T_{r), namely 722 (= 733) and Fi, are illus- 
trated in Fig. 1 as functions of x and ?/ for cr = 5, i? = 1 and v'^ G {0.3,0.6,0.9}. The corresponding plots of 
722 and Fi versus z arc identical to those versus y. The constitutive parameters 722 and Fi arc continuous 
functions of r for all values of e [Ojl)- Furthermore, the constitutive-parameter space is approximately 
partitioned into two disjoint regions with 722 and Fi being approximately constant-valued in each. That is, 
we have 

(i) an inner region — which corresponds to the warp bubble of Alcubierre spacetime — wherein 722 ~ 
1 — {v'^) and Fi ~ w^, and 

(ii) an outer region wherein 722 ~ 1 and Fi w 0. 

The inner region is shaped like a prolate spheroid whose major axis is aligned parallel to the x axis. The 
prolate spheroid becomes increasingly elongated as the relative speed v'^ increases. 

For the presentation of ray trajectories in Sec. 4, let us introduce the semi-major axis length um and 
semi-minor axis length am of the ellipse representing the inner region in the xy plane, defined via 



x-T{aMx) = -£(0) 



(14) 



3 Analysis of quasi-plane waves 

As a precursor to our investigation of ray trajectories, we first consider a quasi-plane wave whose electric 
and magnetic fields are of the form [24] 



E{r,t) = Rc{E„{r)c^p[i{kok-r-iut)]} 
H{r, t) = Rc { ii„(r) exp [i {k„k ■ r - tot)]} 



(15) 



The quantities E_^[r) and H_^{r) in eqs. (15) arc spatially varying, complex-valued amplitudes; w is the 
angular frequency; and the wavenumber kg = uj^Je^^. Within the quasi-planewave regime the relative 
wavevector k varies with r, but it is convenient to omit the dependency on r in our notational representation 
of k. 

The Maxwell curl postulates in combination with the constitutive relations (8) and electromagnetic fields 
(15) yield the nonhomogeneous vector differential equations 



[V (fc • r) - r(r)] X E,Xr) - .[^lir) ■ H,{r) = -t^ V x E„{r) 
[V {k-r)- Ul)] X H,{r) + .[^l{r) ■ E„{r) = — ^ V x H„(r) 

V Mo — *"'0 



(16) 



Under the geometric-optics approximation, the constitutive parameters are assumed to vary only very slowly 
over the distance of a wavelength. Thus, V (fc • r ) « fc and the right sides of eqs. (16) are approximately 
null-valued. Hence eqs. (16) reduce to [25] 

^etji_^x,y)-p-j{^x,y)-p L + PP- j{.x,y)^ ■ E^{r) ^ (17) 

wherein the vector p ^ k — r(a;, y) is introduced. The existence of nonzero solutions to eq. (17) imposes the 
condition 

'H = detj{x,y)-p-'y{x,y)-p = 0. (18) 
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Equation (18) represents the dispersion relation from which the magnitude k of the relative wavevector 
k may be extracted as follows. Writing k = k {sin 6 cos(f> x + sm9 sin (j)y + cos6 z), we see that the left side 
of eq. (18) is quadratic in k; hence, 



k=k^( 



-b± \J\P' - iac 
2^ ' 



where the coefficients 



1 - {lif] (cos2 e + sin^ esin^ 0) + {l - [1 - f{r,)f (v'f] 



sin^ 6 cos^ I 



first {v',f 



2k, 



b = -- 



Wsfirs) {l + [1 - /(r«)] «)'} sine cos 



first {v',f 



k?Al- 



[l + firs)[2-firs)iv'sf]]iv',f) 



l-[l-firstiv'f 



(19) 

(20) 
(21) 
(22) 



Let us note that fc+(7r — 0, tt + </>) = — A;~ (0, 0); i.e., k~ is the wavenumber for a quasi-plane wave travelling 
in the opposite direction to the quasi- plane wave with wavenumber fc"'" . This observation is a manifestation 
of the unirefringence of vacuum [26] . 

The discriminant term b^ — Aac in eq. (19) can have a negative value. Therefore the relative wavenumber 
k may be complex-valued with nonzero imaginary part, in spite of the fact that all the coefficients of the 
dispersion relation (18) have real values. However, Im{fc} 7^ is indicative of evanescent waves which we 
exclude from our study of ray trajectories. 

The partition of the fc-phase space into a propagating-wave regime and an evanescent-wave regime is 
illustrated in Fig. 2, wherein the directions of k for which Im {fc} ^ arc represented at the coordinate origin 
for v'g S {0.60,0.62,0.70,0.90}. We choose the parameter values (7 = 5 and i? = 1 for these plots. Only 
the directions in one octant of the unit sphere need be displayed, because of symmetry. At relative speeds 
Wg < 0.6, the relative wavenumbers are wholly real- valued for all propagation directions. As v'^ increases just 
beyond 0.6, relative wavenumbers with nonzero imaginary parts emerge for k directed in the yz plane. As 
v'g increases further, Imjfc} ^ occurs at increasingly larger values of x • fc; in the limit v'^ ^ 1, we find that 
Im {k} 7^ occurs for all directions of propagation. The same trend is observed at locations throughout the 
inner region of the constitutive-parameter space referred to in our discussion of Fig. 1. In the outer region, 
however, k is everywhere real- valued for all propagation directions. 



4 Ray trajectories 



A convenient Hamiltonian function for our ray-tracing study is provided by the scalar quantity % introduced 
in eq. (18). We parameterize the ray trajectories in terms of r via ^(t); similarly, the parameterization fc(r) 
is used for the relative wavevector. The ray trajectories are thus governed by the coupled vector differential 
equations [27, 28] 

rlr 



dr 
dr 

wherein the shorthand Vg = xd/dqx + yd/dqy 

appropriate initial conditions r(0) and fc(0) have been specified, the system (23) can be solved for r(T 
kir) using standard numerical methods — e.g., the Runge-Kutta method [28]. 

That the direction of a ray trajectory, as provided by the direction of V^H, is aligned with the direction 
of energy flux (quantified as the time-averaged Poynting vector) was demonstrated previously for a general 



for q = qx X + q.y y + qz z_ is employed. 



(23) 



Once 
) and 
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Tamm medium, under the proviso that 7(r) is either positive-definite or negative-definite [25, 29]. It is 

clear from eq. (10) that botli distinct eigenvalues of j{r_), namely 1 and yi , are positive-valued. 

Hence, the direction of energy flux is indeed aligned with the ray trajectories deduced from eqs. (23). 

We begin our presentation of ray trajectories by considering an array of rays in the xy plane, initially 
parallel to the x axis and equally spaced. That is, we take r(0) = xoX + yoy with xq fixed and —1.5am < 
2/0 < 1-5 «„. As in Figs. 1 and 2, we set cr = 5 and R= 1. In Fig. 3, ray trajectories are shown for xq < 0, 
fc(0) = X, and v'^ £ {0.3, 0.6, 0.9}. The inner region of the constitutive-parameter space referred to in our 
discussion of Fig. 1 is shown as a shaded (yellow) ellipse (with semi-major axis length um and semi-minor 
axis length Um, per eqs. (14)), which becomes more eccentric at larger values of v'^. The inner region is seen 
to have a focusing effect, with the focus lying on the +x axis. Furthermore, the focus shifts towards the 
coordinate origin as the relative speed increases. 

That the Tamm medium is not reciprocal in the Lorentz sense [30] is vividly illustrated by a comparison 
of Figs. 3 and 4. The scenario represented in Fig. 4 is the same as that of Fig. 3 except that Xq > and 
A;(0) = —X. Quite unlike Fig. 3, there is no evidence of focusing by the inner region in Fig. 4. On the 
contrary, rays appear to diverge as they pass through the inner region at low values of v'^, while rays are 
almost entirely excluded from the inner region altogether at v'^ = 0.9. 

Further insight into the absence of Lorentz-reciprocity of the Tamm medium is provided in Fig. 5 wherein 
ray trajectories initially parallel to the y axis and equally spaced are presented. The initial relative wavevector 
for these trajectories is fc(0) = y, and we have set r (0) = xq x + yoy with fixed yo < while —1.5 um < xq < 
1.5 um- As for Figs. 1-4, a = 5 and R = 1. The plots in Fig. 5 are clearly asymmetric with respect to the y 
axis, and the ray trajectories become progressively excluded from the inner region as v'^ increases. 

Rays initially propagating in radial directions in the xy plane are represented in Fig. 6 for cr = 10 and 
R ~ 1. We track rays which emanate from point sources in the inner region (at the coordinate origin) 
and in the outer region at locations on the positive and negative x axis; i.e., r(0) = xqx. The relative 
speed Wg S {0.3, 0.6, 0.9}. Equally-spaced angular directions for the initial relative wavevector fc(0) were 
considered. However, some initial directions in the inner region correspond to evanescent waves and these 
are not represented in Fig. 6. The proportion of directions which correspond to evanescent waves increases 
as Vg increases for sources in the inner region. Indeed, for v'^ = 0.9 with the source at the coordinate origin, 
only 30% of the possible fe(0) radial directions correspond to propagating rays. The general trends apparent 
in Figs. 3-5 are also apparent in Fig. 6. That is, the inner region has a focusing effect for sources located 
outside the inner region with xq < 0; for sources located outside the inner region with xq > ray trajectories 
tend to be progressively excluded from the inner region as v'^ increases. 

Let us now turn to the influence of the thickness of the transition zone between the inner and outer 
regions, as dictated by the parameter a via the scalar function /. In Fig. 7 ray trajectories are provided 
which correspond to the scenario of Fig. 6 but with a ~ 1 and 25. We note that a ~ 25 results in a more 
sharply defined top-hat profile with straightcr sides for /, whereas a = 1 results in a profile more rounded 
sides, as compared to cr = 10 which was used for Fig. 6. Comparing Figs. 6 and 7, we deduce that, although 
the change in the direction of rays at the boundary between the inner and outer regions becomes more 
pronounced as a increases, the general pattern of ray trajectories remains largely unaffected. 

For clarity of representation, ray trajectories restricted to the xy plane were considered in Figs. 3-7. 
Trajectories of the same form can be observed in the xz plane. The trajectories for the yz plane are likewise 
similar, albeit then the inner region of the constitutive-parameter space is obviously circular in shape, 
regardless of the relative speed v'^. 

5 Closing remarks 

A flat-spacetime representation of the Alcubierre spacetime has been established by means of a Tamm 
medium which is asymptotically identical to vacuum and has constitutive parameters which are continuous 
functions of the spatial coordinates. Thus, the Tamm medium is amenable to physical realization as a 
nanostructurcd mctamatcrial. An alternative approach — which utilizes a Galilean transformation instead 
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of the Lorentz transformation (4) — gives rise to a Tamm medium which is not asymptotically identical to 
vacuum and is accordingly less well-suited to physical realization [31]. 

Our geometric-optics study has revealed that ray trajectories are acutely sensitive to the relative speed 
v',. of the warp bubble, and to the direction of its velocity i;^, but rather less sensitive to the thickness of 
the transition zone between the warp bubble and its background. In particular, for rays which travel in the 
same direction as x/g, the warp bubble acts as a focusing lens, especially at large values of v'^. 
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support of his research activities. 
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Figure 1: The constitutive parameters 722 and Fi plotted versus x for y = z = 0, and versus y for x = z = 0. 
Parameter values: ct = 5, i? = 1 and v'^ — 0.9 (blue, broken dashed curves), 0.6 (green, dashed curves) and 
0.3 (red, solid curves). 
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Figure 2: Maps illustrating the directions oi k ~ k/ k = kx x + ky ij + kz z_ tor which Im {k} = (blue) and 
Im{fc} 7^ (red), at the coordinate origin for v'^ G {0.60,0.62,0.70,0.90}. Parameter values: (7 = 5 and 
R=l. 
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v's=0.3, xo=-28.0 v's=0.6, xo=-28.0 
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Figure 3: Trajectories for rays in the xy plane, initially parallel to the x axis and equally spaced. That 
is, r(0) = Xox + yoy with fixed xq = —28 and — 1.5a„i < y^ < 1.5 a™, fc(0) = x and v'^ G {0.3,0.6,0.9}. 
Parameter values: cr = 5 and i? = 1. 
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Figure 5: As Fig. 3 except that r(0) = xox + yoy with — 1.5aA/ < xq < 1.5 om and fixed yo — —28, and 
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Figure 6: Trajectories for rays in the xy plane, emanating from sourees on the x axis, i.e., r(0) = x^x^ at 
equally spaced angular directions of fc(0), for v'^ S {0.3,0.6,0.9}. Parameter values: cr = 10 and R = l. 
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